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Abstract 


Active Galactic Nuclei (AGN) may emit highly collimated and intense jets 
of relativistic electrons which upscatter ambient photons. These electrons can 
also scatter off the cold dark matter halo of the galaxy to produce high energy 
photons which have a more isotropic signature than the up-scattered photons 
from QED processes. We propose to look for these high energy photons com¬ 
ing from AGN as a method to detect dark matter. As a primary example we 
work out the expected signal from electrons scattering off the lightest super- 
symmetric partner into a photon plus selectron. Using the optimistic side of 
astrophysical uncertainties, we still find the signal from M87 or Gentaurus A, 
two close-by AGN, smaller than the sensitivities expected of the currently pro¬ 
posed photon detectors. However, long running photon detectors and future 
detectors of higher sensitivity might be able to distinguish a signal from AGN 
sources. In order to have confidence that new physics sources are discernible, 
we also emphasize the importance of multi-wavelength studies of AGN with 
varying jet axis orientation to the earth. 
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Astrophysics observations and cosmological argnments point to the necessity for a snb- 


stantial amonnt of cold dark matter in the nniverse [j^. It is even possible that the invisible 
cold dark matter constitntes nearly critical density (pc = Sif^/SvrG) in the standard = 1 
big-bang cosmology. Dark matter detection to date has only been accomplished by witness¬ 
ing gravitational effects (rotation cnrves of galaxies, Hnbble flow distortions, gravitational 
lensing, etc.). Since the dark matter is probably not electrically charged or color charged [@, 
its most likely non-gravitational interactions are dne to the weak force. Hence, the hypothe¬ 
sis of stable weakly interacting massive particles (WIMPs) has been considered [Q]. We will 
nse the bino (snperpartner of the hypercharge gange boson) as the lightest snpersymmetric 
partner (LSP) in onr example scattering process. 

Experimental effort is now nnderway on many fronts |p to see dark matter via its weak 
interactions. Cryogenic table top experiments which measnre nnclear recoil of ambient LSPs 
interacting with Silicon, Germaninm, and other elements are progressing. One can also 
search for nentrino yields from the annihilation of captnred LSPs in the core of the snn or 
earth. And searches are underway to look for positrons, electrons, protons and antiprotons 
from the annihilations of LSPs in the galactic halo. 

Here, we propose another idea to search for LSPs: photons originating in the hnal state 
of electron-LSP scattering near the center of Active Galactic Nuclei (AGN). The principle 
property of an AGN of the Blazar type which makes this idea potentially feasible is the 
high energy and high flux electron beams which are thought to be emanating from their 
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centers The high energy jets are thought to arise from a highly relativistic shock wave 
traveling along the perpendicular axis of an accretion disk associated with a massive black 
hole at the center (10® — 10® Mq). Shock waves generated by the central engine accelerate 
electrons and positrons which in turn produce high energy photons mainly through inverse 
Compton scattering of the UV radiation from a thermal distribution of photons originating 
external to the electron/positron jet in the accretion disk |^. Other processes that could be 
important include inverse Compton scattering of synchrotron radiation, and pair annihilation 
of the electrons and positrons into high energy photons. We will loosely call all these 
processes “up-scattering” of high energy electrons and positrons into high energy photons. 

The lab frame, which is the accretion disk frame, or equivalently our observation frame 
on earth (for close-by AGN), has a distribution of electrons which scales roughly as E~‘^, 
and is highly collimated with a beaming factor (solid angle fraction) of about 10“^ |^. Since 
the upscattered radiation can be measured, the coefficient out in front of the assumed E~‘^ 
can be measured to obtain an estimate of the electron/positron flux ejected from the AGN. 

The total power of the AGN in photon radiation is approximately AE^ = 10“^^/ — 
10^®/ ergs/s in the energy range of 100 MeV < E-y < 5 GeV 0, where / is the beaming 
factor (fraction of solid angle) of the electron beam. We shall take / ~ 10“^ in accordance 
with [^. Since the hadronic component of the AGN jet is probably less important Q to the 
photon radiation than the electron/positron constituents, we will make the approximation 
that all the photon radiation power is converted from the electron jet (often we use “electron” 
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to mean both electrons and positrons). This enables us to estimate the intensity and energy 


distribution of the electron beam. We hnd the approximate value by using the same 
energy prohle for the electrons as the photons: 


dEe 



( 1 ) 


We can estimate the value of no to be 


no ~ 


1 

ml logs GeV/lOOMeV 


~ lO^V GeV/s. 


( 2 ) 


Given this normalization, and given this distribution between 1 MeV to 10 TeV, the total 
amount of energy ejected from the AGN is about one solar mass per year. It might be that 
much more energy is ejected from the AGN, and that all the electron beam is not transmitted 
into photons, in which case the above estimate for no would be too low. However, it should 
also be kept in mind that there are other models of high energy photon jets from AGNs 
which do not require such intense beams of electrons 0. 

If all the electron power went into merely up-scattering photons through simple QED 
processes, then there would be no electrons left to interact with dark matter. Below it will 
become apparent that the penetration depth of the electrons is not crucial unless it is less 
than a few parsecs from the central engine, since almost all the electrons have to interact with 
the LSPs within about ten parsecs of the accretion disk in order for a signal to be viable at 
all. Together with the observation that a relatively small fraction of the electrons do scatter 
off the LSPs, we are free to use Eq. as our particle beam source which interacts with the 
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dark matter. Although there are some reasons to believe that the electrons could upscatter 


the photons far away from the central engine |^, the reader should keep in mind that a 
detailed model description of the electron/photon beam will involve a distance dependent 
attenuation factor for the electron beam. Current models are not precise enough to calculate 
the attenuation so we have assumed optimistically that it is negligible. 

The dark matter target distribution must also be modeled [^,0. Several differing pro¬ 
posals are in the literature. Here, we choose to model the dark matter distribution after 


Ref. which postulates that p ~ r We choose this model over the more traditional 


distributions of p ~ l/(r^-|-r^), where Vc is a hard-core density radius which smoothes p near 
r = 0, for two reasons: (1) The lack of a hard-core density radius in this model is becoming 


more and more observed in numerical simulations , and (2) large enhancements of p near 


r ~ 0 are needed for our proposed signal to be interesting. That is, without a cusp in the 


dark matter distribution [O, or equivalently, without Vc < few parsecs, our resulting signal 


flux of photons would be too low to be detected. 


To be explicit, we use the dark matter distribution 


p(r) = Pa - 
\r 


1.8 


( 3 ) 


where a is some arbitrary distance from the center of the galaxy and pa is the local density 
at that distance. For numerical purposes we will use a = 8 kpc and pa = 0.3GeV/cm^, 
in accord with the density profile of our galaxy. It is important to estimate the minimum 


radius size, at which this density function breaks down. The different physics sources 
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central black hole would yield the largest of ~ Ipc (for M ~ 10®Mq). The assumption 


of a black hole at the center of the galaxy is quite applicable to our AGN study, and we 
adopt the black-hole hypothesis as the origin of The capture radius of dark matter is 
then between about 0.01 pc to 10 pc depending on black hole mass, dark matter density, 
etc. For numerical purposes we will use both 1 pc and 1 kpc. 

We now have modeled both the target (dark matter distribution) and the source (high 
energy electron beam), and so we can estimate the flux of photons originating from e^x 
e ^7 scattering. Our notation identihes the bino with x and the right-handed selectron 
(superpartner of the right-handed electron) as cr. Feynman diagrams for the process are 
given in Fig. |^. 

Before setting up the differential flux, it is instructive to summarize some basic relativistic 
kinematics and notation. For close-by AGN, the lab frame of the electron/LSP collision is 
the observer frame. Thus E^. is the incident electron energy, making the center of mass 
energy equal to 


s = 


+ 2.Fe), 


(4) 


and the resulting photon energy in the lab frame is 


s — 


( 5 ) 


E. 


'7 


2m^ 


In the c.m. frame (primed quantities) the electron and photon energies are 
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( 6 ) 


^ s-ml ^ s-ml^ 

^ 2 ^ ^ 2yi • 

Relativistic transformations of fonr vectors in the lab frame to the bino/electron c.m. frame 
are carried ont by 


P = 




(7) 


Ee + 

The angle x = cos 9 with respect to ;§ (jet axis) in the lab frame is related to the angle 


x' = cos 6' in the c.m. according to 


X = 


P + x' 

1 + Px'' 


( 8 ) 


In the c.m. frame the differential cross-section for ex ey scattering in the limit that 


s S> rriy, rriej^ simplifies to 


da 


a 


^ 1 /1 + cos 9'' 


dVt' 2 cos^ 9w s \ 1 — cos 9 


/ ; ’ 


(9) 


where a = 1/137 is the QED conpling constant, and 9w is the weak mixing angle. (Dne 
to hnite electron mass and selectron width there is no singnlarity in the cross-section as 
9' —>■ 0.) We keep the masses of y and cr in onr nnmerical work. To cast this in the lab 
frame we nse Eq. | to snbstitnte x for x' and to calcnlate dx'/dx as a fnnction of x. 

The total differential e^x scattering flnx observed at the earth can be parameterized as 


dF ^ events 


dEf, crrPsecGcY d^Qj^ 


da 

dVt 


XO 


dC 

dEP 


( 10 ) 


where xq = cos^o is the fixed angle between the jet axis and the axis which points from 
the AGN to the earth. Also, d^cN is the distance to the AGN and dC/dEf. is the effective 


Inminosity given by 
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( 11 ) 


dC ^ df (pj^) 
dEg dEe 


The factor (p^^x) average of the r dependent dark matter density function times the 

effective length that the electrons pass through, and is dehned as 


I * ! t HJjtjU 

{px^x) = / drp{r) ~ r^i^pa 




( 12 ) 


where the latter approximation is for rmax S> Tmin- 


There is a one-to-one correspondence between the electron energy Eg and the photon 


energy E^ given by 


Ee = 


nil — m? + 2myE. 




2m. 


(13) 


Since dE^/dE^ = 1 the differential photon flux is given simply by substituting Ef. of Eq. 


into Eq. IM The integral photon flux is then 


E{Ey) = # photons cm ^sec ^ = / 

J E 


OO (jp 

dE'i 


^dE'X 


(14) 


Figure 2 shows the results for the AGN M87 which is approximately 12 Mpc away, and 


has a jet axis oriented at 42 ± 5 degrees to us [13,13. We also show the result for Centaurus 


A, which is ~ 2.5 Mpc away and is oriented at ~ 68 degrees |^. The supersymmetric 


model used to construct Fig. 2 is a pure bino LSP with mass 60 GeV, and a right-handed 
selectron with mass 100 GeV. The upper solid line is with = Ipc and the bottom solid 
line is with = 1 kpc. The top dashed line is the the 5a sensitivity of the proposed space 


based gamma ray telescope GLAST [Q , and the lower right line is the expected sensitivity 
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of the proposed ground based gamma ray telescope array VERITAS [17]. For GLAST, the 


sensitivity corresponds to one year scanning mode in which the instrument is always pointed 
outward from the earth (approximately 30% duty cycle). For VERITAS, the sensitivity is 
for one week observing time on the source. As we can see, the signal is lower than the 
present sensitivities of these proposed photon detectors. 

The AGN beam as modeled in has a non-zero gamma ray flux at large angles from the 
jet axis. For Centaurus A the flux at 68 degrees off jet axis (i.e., earth direction is 68 degrees 
from jet axis) is expected to be negligible; the electron/positron beam is highly collimated, 
dramatically decreasing the overall flux normalization and especially decreasing the gamma 
ray energy cut-off from ordinary QED processes with higher observing angle. However, for 
M87 whose jet axis is only 42 ± 5 degrees with respect to the earth direction, the intrinsic 
background is substantial and swamps the expected flux from electron/bino scattering up 
to approximately 150 GeV. In the hgure we have plotted this background estimate using 
the same photon beam flux used earlier in deriving the dark matter signal, and we have 
also assumed a 45 degree jet axis orientation for easier comparison with ref. [^. Measuring 
gamma ray fluxes from intermediate angle (10° ^ 6 ^ 50° off jet axis) AGN such as M87 are 
still useful because they can increase our understanding of AGN jets, and possibly the model 
of ref. in particular. Multi-wavelength measurements of the photon flux from a number 
of AGN with different jet axis orientations is highly desirable to more fully characterize the 


AGN beams including determination of the angle of the jet axis to the earth using radio 
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measurements. Such measurements will allow deeper understanding of AGN beams and 


possibly allow us to have conhdence that new physics is discernible. 

Although we typically chose parameters such as r^in and no which gave maximal contri¬ 
bution to the signal, there are many astrophysical uncertainties involved in the calculation 
that could make the signal perhaps larger or smaller than we have estimated here. Summing 
over all large angle Blazar sources will of course help increase the signal over the diffuse back¬ 
ground. Such techniques have already been shown to have a signihcant impact on sensitivity 


in the X-ray region |18 . 


We have focused on the case of a pure bino dark matter candidate. Certainly other 
particles are around, even within supersymmetry, that yield good dark matter candidates. 
These candidates could of course have much different scattering amplitudes than the bino. 


For example, the higgsino dark matter candidates discussed in |T^ would not allow a direct 


coupling of the higgsino to an electron and selectron. One useful hnal state is a chargino and 
a neutrino produced from t channel W exchange. The decay chain from the chargino then 


can include a photon from X 2 Xi7- Such decays could lead to structure in the photon 


energy spectrum, providing an additional signature. In any event, the hnal cross-section 
will have electroweak strength and should not be signihcantly different than the bino case 


studied above. 


We have found that our estimate of the signal for electron-bino scattering into photons 


from a single AGN source is too low to be discerned by currently proposed experiments. 
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but future detectors of much higher sensitivity might be able to see a signal. The limit on 
sensitivity to the photon signal of the most sensitive of the currently proposed experiments 
- ground based Cherenkov telescope arrays - is largely due to the difficulty of subtracting 
pseudo-backgrounds primarily from electron showers in the atmosphere. Therefore, if a 
signal does exist it might be possible to reach the sensitivities required by the continuing 
progress to eliminate these backgrounds, and dramatically extend the observing time on 
candidate sources. Obviously this presents great challenges. 
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FIGURES 



FIG. 1. Feynman diagrams contributing to e^x l^R scattering. 
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FIG. 2. Integral flux of photon signal from the Gentaurus A and M87 Active Galactic Nuclei. 
The top solid (dash-dot) line corresponds to = 1 pc in the dark matter density profile for 
Gentaurus A (M87), and the bottom line corresponds to = 1 kpc. The dashed lines represent 
the sensitivity of the different photon detectors GLAST (one year scanning mode) and VERITAS 
(one week on source). The dotted line represents an estimate of the background flux from the M87 
jet. Gamma ray jet background from Gentaurus A is expected to be negligible. See the text for 
further discussion. 
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